o 



Oh 
I 

o 



> 

o 

o 
o 



% 



Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 27 July 2010 (MN I*TeX style file v2.2) 

Magnetospheric accretion and spin-down of the 
prototypical classical T Tauri star AA Tau 

J.-F. Donati^*, M.B. Skelly\ J. Bouvier^, S.G. Gregory^ K.N. Cranking 
M.M. Jardine^ G.A.J. Hussain^ F. Menard^, C. Dougados^, Y. Unmh^ 
S. Mohanty^, M. Auriere^, J. Morin^'^, R. Fares^ & the MaPP collaboration 

1 LATT-UMR 5572, CNRS & Univ. de Toulouse, 14 Av. E. Belin, F-31400 Toulouse, France 

2 LAOG-UMR 5571, CNRS & Univ. J. Fourier, 414 rue de la Piscine, F-38041 Grenoble, France 
^ School of Physics, Univ. of Exeter, Stocker Road, Exeter EX4 4QL, UK 

^ Crimean Astrophysical Observatory, Nauchny, Crimea 3344^3, Ukraine 

' School of Physics and Astronomy, Univ. of St Andrews, St Andrews, Scotland KY16 9SS, UK 

' ESQ, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 

Department of Physics, Imperial College London, London SW7 2AZ, UK 
' Dublin Institute for Advanced Studies, School of Cosmic Physics, 31 Fitzwilliam Place, Dublin 2, Ireland 



2010 June, MNRAS, submitted 



ABSTRACT 

From observations collected with the ESPaDOnS spectropolarimeter at the 
Canada-France-Hawaii Telescope (CFHT) and with the NARVAL spectropolarimeter 
at the Telescope Bernard Lyot (TBL), we report the detection of Zeeman signatures 
on the prototypical classical T Tauri star (cTTS) AA Tau, both in photospheric lines 
and accretion-powered emission lines. Using time series of unpolarized and circularly 
polarized spectra, we reconstruct at two epochs maps of the magnetic field, surface 
brightness and accretion-powered emission of AA Tau. We find that AA Tau hosts 
a 2 — 3 kG magnetic dipole tilted at ~ 20° to the rotation axis, and of presumably 
dynamo origin. We also show that the magnetic poles of AA Tau host large cool spots 
at photospheric level and accretion regions at chromospheric level. 

The accretion rate at the surface AA Tau at the time of our observations (esti- 
mated from the emission in the He i D3 line mainly) is strongly variable, ranging from 
—9.6 to —8.5 and equal to —9.2 in average (in logarithmic scale and in Mq yr~^); this 
is an order of magnitude smaller than the disc accretion rate at which the magnetic 
truncation radius (below which the disc is disrupted by the stellar magnetic field) 
matches the corotation radius (where the Keplerian period equals the stellar rotation 
period) ~ a necessary condition for accretion to occur. It suggests that AA Tau is 
largely in the propeller regime, with most of the accreting material in the inner disc 
regions being expelled outwards and only a small fraction accreted towards the surface 
of the star. The strong variability in the observed surface mass-accretion rate and the 
systematic time-lag of optical occultations (by the warped accretion disc) with respect 
to magnetic and accretion-powered emission maxima also support this conclusion. 

Our results imply that AA Tau is being actively spun-down by the star-disc mag- 
netic coupling and appears as an ideal laboratory for studying angular momentum 
losses of forming Suns in the propeller regime. 

Key words: stars: magnetic fields - stars: formation - stars: imaging - stars: rotation 
- stars: individual: AA Tau - techniques: spectropolarimetry 
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1 INTRODUCTION 

Classical T Tauri stars (cTTSs) are young low-mass stars 
still contracting towards the main sequence and surrounded 
by gaseous and dusty accretion discs; they represent the im- 
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portant formation stage that stars with masses lower than 
~ 3 M© undergo at an age of a few Myrs, and during which 
they build up their exotic planetary systems. 

Observations reveal that magnetic fields play a crucial 
role at the cTTS stage. Thanks to strong large-scale fields 
(e.g., iDonati fc Landstreetl 12003 . for a review), cTTSs are 
able to evacuate the core regions of their accretion discs, 
and to connect to the inner rim of their discs via discrete 
accretion funnels or veils through which mate rial is accreted 
and a ngular momentum is dissipated (e.g., iBouvier et al.l 
l2007al . for a review). Accretion discs are also expected to 
be magnetic, with fields enhancing accretion rates, generat- 
ing powerful jets and modifying planet formation/migration 
mechanisms. Surveys of the magnetic properties of cTTSs 
and their accretion discs are thus critically needed to under- 
stand the early history of low-mass stars in general and of 
the Sun in particular. 

Magnetic Protostars and Planets (MaPP) is an inter- 
national project focussing specifically on this issue. It has 
been granted 690 hr of observing time over 9 consecu- 
tive semesters (2008b to 2012b) with the ESPaDOnS spec- 
tropolarimeter on the 3.8-m Canada-France-Hawaii Tele- 
scope (CFHT) to survey 15 cTTSs and 3 protostellar ac- 
cretion discs of FU Ori type (FUOrs); it also regularly 
benefits from contemporaneous observations with the NAR- 
VAL spectropolarimeter on the 2-m Telescope Bernard Lyot 
(TBL) as well as photometric observations from Crimea, 
Uzbekistan and Armenia. Additional multiwavelength ob- 
servations from space (e.g., XMM-Newton/Chandra/HST) 
and/or from the ground (e.g., HARPS) are also organised in 
conjunction with MaPP campaigns on a few specific targets, 
providing deeper insights into the physical processes under 
scrutiny (and in particular magnetospheric accretion). 

MaPP primarily collects spectropolarimetric data prob- 
ing the large-scale magnetic fields of cTTSs through the 
Zeeman signatures they generate in photospheric line pro- 
files; it also allows the detection of Zeeman signatures in 
the accretion spots located at the footpoints of accretion 
funnels. By monitoring these Zeeman signatures over sev- 
eral successive rotation cycles (to filter out intrinsic vari- 
ability and to retrieve rotational modulation more reliably 
and efficiently), MaPP can reconstruct maps of the large- 
scale magnetic fields of cTTSs and simultaneously recover 
the location of accretion spots. By extrapolating from sur- 
face magnetograms, one can finally obtain an approximate 
3D description of the magnetosphere, allowing more realistic 
and quantitative models of the geometry of accretion funnels 
and more generally of the star/disc magnetic coupling. Ini- 
tial pre-MaPP stu dies of a few cTTSs were presented in sev- 
eral publications (IDonati et al.l 20071. 2008a : IJardine et al 



20081 : iGregorv et al.ll2008l : iHussain et al.ll2009l : lDonati et al 



20101 1 to validate the main assumptions underlying the imag- 



ing code and demonstrate the overall feasibility of the mod- 
elling. 

The present paper is the first in a series dedicated to 
MaPP data and results, using a more general and better 
suited set of modelling assumptions and a more mature ver- 
sion of the imaging code; ultimately, the main goal is to pro- 
vide estimates and statistics on the topologies of large-scale 
magnetic fields of cTTSs and on the locations of their accre- 
tion spots, allowing the diagnosis of how they correlate with, 
e.g., stellar masses and rotation rates (in a way similar to 



that achieved on main-sequence stars, IDonati fc Landstreeti 
[2OO9I), but also on parameters more specific to cTTSs, e.g., 
ages and accretion/ejection properties. This first MaPP 
study concentrates on the prototypical cTTS AA Tau (see 
Sec. [2] for a quick summary of the main stellar parameters 
relevant to this study). We report here spectropolarimetric 
and photometric observations of AA Tau (Sec. [3} , describe 
the variations of photospheric lines and accretion proxies 
(Sec. 131) and their subsequent modelling (Sec.[S|. We finally 
discuss the implications of these new results for our under- 
standing of magnetospheric accretion processes in cTTSs 
(SecO. 



2 AA TAU 

AA Tau is a well-known cTTS showing strong Ha emission 
and IR excesses demonstrating the presence of a gaseous 
and dusty acc retion disc surrounding the protostar (e.g., 
iBouvier et "all [1999). A fit to the observed BVRI photo- 
metric fluxes (showing B — V, V — Re and V — Ic col- 
ors of 1.5, 0.9 and 1.9 respectively at times of maximum 
brightness and min imum accretion, see Figs. 3 and 6 of 
IBouvier et al.l |2003| ) suggests that AA Tau has a photo- 
spheric temperature of Teg ~ 4000 K an d a visual redden- 
ing of Ay ~ 0.8 mag (|Bouvier et al.lll999l), i mplving a visual 
bolometric correction of —1.8 i Bessell et al.lll998. 1. Given its 
distance (~ 140 pc) and maximum vis ual brightness (corre - 
sponding to a magnitude of ~12.3, e.g.. lGrankin et al.ll2007l ). 
AA Tau has a bolometric luminosity equal to that of the 
Sun (within about 0.1 dex) and thus a radius of 7?* ~ 2 Rq, 
suggestin g a mass of Mj, ~ 0.7 M© and an age of about 
1.5 Myr llSiess et al.ll2000l) in ag reement with previous esti- 
mates (e.g.. IBouvier et al.lll999r) . This implies in particular 
that AA Tau is fully convective and very similar (in mass, 
radius, a ge and rotation rate) to the other prototypical cTTS 
BP Tau (JDonati et aLll2008al '). 

AA Tau is known to undergo periodic eclipses with ir- 
regular depths, likely caused by partial occul tations by a 
warp e d accretion d isc viewed close to edge-on (jBouvier et al.l 
1 19991 . I2OO3I . l2007bl ). As the warp in the accretion disc is 
presumably caused by the magnetic field of the star, the 
recurrence period of these eclipses, equal to 8.22 ± 0.03 d 
(JBouvier et al.ll2007bl ). is actually tracing the rotation pe- 
riod of the star. Given that _R* ~ 2 Rq, the estimated ro- 
tation period implies that the equatorial rotation velocity 
of AA Tau is 12.3 kms~^; the measured line-of-sight pro- 
jected rotation velocity (denoted vsini where i is the incli- 
nation of the rotat ion axis to the line of sight, and equal 
to 11.3±0.7 kms~^ IBouvier et al.ll2003l ) independently con- 
firms that AA Tau is viewed from the Earth at an inclination 
angle of about 70 ± 10°. With these parameters, the radius 
at which the Keplerian period is equal to the rotation period 
of the star, called the corotation radius and denoted rcor, is 
equal to ~ 7.6 i?* or 0.07 AU. 

Mass accretion on AA Tau is reported to be smaller 
than the average rate expected for cTTSs of similar masses 
(e.g., Johns-KruU 2007). The spectrum of AA Tau exhibits 
all the usual accretion proxies, in particular Ha, H/3, He I 
-D3 and Ca 11 infrared triplet (IRT) emission. From their 
strength (in particular that of He I emission) and the cor- 
responding line fluxes (see Sec. |3]l in our spectra, and using 
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empirical correla tions from the published literature (e.g., 
iFang et al.l 120091 ) . we can estimate the logarithmic mass- 
accretion rate at the surface of A A Tau (in M© yr~^), found 
to vary from —9.6 to —8.5 throughout our runs (e.g., in 2008 
December) and equal to about —9.2 in average. This is up to 
10 times smaller than the accretion rate of BP Tau (equal to 
—8.6 when estimated with the same method, iDonati et al.l 
[2010!). Accretion on A A Tau is known to be intrinsically 
variable on a time scale of a few days (e.g., iBouvier et al.l 
I2OO3I . l2007bl ). Optical veiling, i.e., the apparent weakening 
of the photospheric spectrum (presumably caused by accre- 
tion), is also often ob served at a modera t e (and variable) 
level on AA Tau (e.g.. lBouvier et al.|[20o3 . l2007bl 'l. 



3 OBSERVATIONS 

Spectropolarimetric observations of AA Tau were collected 
in 2008 December and 2009 January using ESPaDOnS on 
the CFHT. ESPaDOnS collects stellar spectra spanning the 
whole optical domain (from 370 to 1,000 nm) at a resolving 
power of 65,000 (i.e., 4.6 kms^^) and with a spectral sam- 

J iling of 2.6 k ms~^, in either circular or linear polarisation 
Donay I2OO3) . A total of 18 circular polarisation spectra 
were collected in 2 separate blocks shifted by about 1 month, 
with 11 spectra over a period of 15 d in 2008 December and 
7 spectra over a period of 8 d in 2009 January; all polari- 
sation spectra consist of 4 individual subexposures lasting 
each 1200 s and taken in different polarimeter configurations 
to allow the removal of all spurious polarisation signatures 
at first order. Six additional spectra were collected a year 
before (in 2007 December and 2008 January) over a period 
of 12 d, with the ESPaDOnS twin NARVAL on the TBL, 
with slightly shorter exposure times. 

All raw frames are processed with Libre ESpRIT, 
a fully automatic reduction package/pipeline available at 
CFHT and TBL. It automatically performs optimal extrac- 
tion of ESPaDOnS unpolarized (Stokes 7) and circularly 
polarized (Stokes V) spectra g r ossly following the proce- 
dure described in iDonati et al.l (|l997h . The velocity step 
corresponding to CCD pixels is about 2.6 kms~^; however, 
thanks to the fact that the spectrograph slit is tilted with 
respect to the CCD lines, spectra corresponding to different 
CCD columns across each order feature a different pixel sam- 
pling. Libre ESpRIT uses this opportunity to carry out op- 
timal extraction of each spectrum on a sampling grid denser 
than the original CCD sampling, with a spectral velocity 
step set to about 0.7 CCD pixel (i.e. 1.8 kms~^). All spec- 
tra are automatically corrected of spectral shifts resulting 
from instrumental effects (e.g., mechanical flexures, temper- 
ature or pressure variations) using telluric lines as a refer- 
ence. Though not perfect, this procedure provides spectra 
with a relative rad ial velocity (RV) pre cision of better than 
0.030 kms"^ fe.g.. lDonati et al.ll2008bl ). 

The peak signal-to-noise ratios (S/N, per 2.6 kms~^ 
velocity bin) achieved on the collected spectra (i.e., the se- 
quence of 4 subexposures) range between 100 and 200 for 
ESPaDOnS data (except for one spectrum recorded in poor 
weather conditions) and between 60 and 90 for NARVAL 
data (directly reflecting the smaller collecting area). Rota- 
tional cycles E are computed from heliocentric Julian dates 



Table 1. Journal of ESPaDOnS/CFHT observations collected in 
2008 December and 2009 January. Columns 1 — 4 respectively list 
the UT date, the heliocentric Julian date and UT time (both at 
mid-exposure), and the peak signal to noise ratio (per 2.6 kms^^ 
velocity bin) of each observation (i.e., each sequence of 4 X 1200 s 
subexposures). Column 5 lists the rnis noise level (relative to 
the unpolarized continuum level /c and per 1.8 kms~^ velocity 
bin) in the circular polarization profile produced by Least-Squares 
Deconvolution (LSD), while column 6 indicates the rotational cy- 
cle associated with each exposure (using the ephemeris given by 
Eq.[l]|. 



Date 


HJD 


UT 


S/N 


0"LSD 


Cycle 




(2,454,000+) 


(h:m:s) 




(10-4/c) 


(1+) 


Dec 06 


806.97106 


11:11:06 


180 


2.5 


48.510 


Dec 07 


807.92124 


09:59:24 


200 


2.3 


48.625 


Dec 08 


808.94040 


10:27:02 


160 


3.0 


48.749 


Dec 09 


809.94201 


10:29:24 


130 


3.7 


48.871 


Dec 10 


810.92881 


10:10:27 


100 


5.8 


48.991 


Dec 15 


815.92556 


10:06:02 


190 


2.4 


49.599 


Dec 16 


816.90946 


09:42:55 


60 


8.7 


49.719 


Dec 17 


817.90170 


09:31:49 


140 


3.7 


49.840 


Dec 18 


818.89873 


09:27:36 


160 


3.2 


49.961 


Dec 19 


819.89824 


09:26:57 


110 


4.6 


50.082 


Dec 20 


820.89756 


09:26:02 


120 


4.3 


50.204 


Jan 07 


838.78706 


06:48:28 


160 


2.8 


52.380 


Jan 09 


840.75479 


06:02:13 


180 


2.6 


52.620 


Jan 10 


841.74685 


05:50:53 


190 


2.3 


52.740 


Jan 11 


842.89978 


09:31:14 


140 


3.8 


52.881 


Jan 12 


843.79311 


06:57:43 


110 


4.9 


52.989 


Jan 13 


844.73006 


05:27:02 


130 


3.9 


53.103 


Jan 14 


845.75784 


06:07:09 


180 


2.7 


53.228 



Table 2. Same as Table [T] for the additional NARVAL/TBL ob- 
servations collected in 2007 December and 2008 January. The 
exposure time of each sequence slightly varies (from 4 X 1000 to 
4 X 1200 s) from night to night. 



Date 


HJD 


UT 


S/N 


0"LSD 


Cycle 




(2,454,000-1-) 


(h:m:s) 




(10-4/c) 


(1+) 


Dec 27 


462.36191 


20:35:13 


70 


10.0 


6.587 


Dec 31 


465.54716 


01:02:15 


80 


7.0 


6.974 


Dec 31 


466.37216 


20:50:19 


90 


6.1 


7.075 


Jan 01 


467.37235 


20:50:41 


80 


7.7 


7.196 


Jan 03 


468.55236 


01:10:01 


60 


8.3 


7.340 


Jan 07 


473.38037 


21:02:49 


70 


7.2 


7.927 



according to the ephemeris: 



HJD = 2454400.0 -f- 8.225 



(1) 



where the rotation period is taken from iBouvier et al.l 
(2007b|). The full journal of observations is presented in Ta- 
bles [T] & [H for ESPaDOnS/CFHT and NARVAL/TBL data 

respectively. 

Least-Squares Deconvolution (LSD. IDonati et al.lll997n 
was applied to all observations. The line list we employed 
for LS D is compute d from an Atlas9 LTE model atmo- 
sphere (|Kurucall993l ) and corresponds to a K7 spectral type 
(Teff = 4,000 K and logg = 3.5) appropriate for AA Tau. 
Only moderate to strong atomic spectral lines (with line- 
to-continuum core depressions larger than 40% prior to all 
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non-thermal broadening) are included in this list; spectral 
regions with strong lines mostly formed outside the photo- 
sphere (e.g., Balmer, He, Ca ll H, K and IRT lines) and/or 
heavily crowded with telluric lines were discarded. Alto- 
gether, more than 9,000 spectral features are used in this 
process, with about 40% of them from Fe I. Expressed in 
units of the unpolarized continuum level 7c, the average 
noise levels of the resulting Stokes V LSD signatures are 
ranging from 2.3 to 5.7xl0~* per 1.8 kms~^ velocity bin for 
ESPaDOnS data (except on Dec 16) and from 6 to 10x10"'' 
for NARVAL data. 

For estimating when eclipses of AA Tau are occurring, 
contemporaneous photometry was collected from Crimean 
Astrophysical Observatory (CrAO) over a period extending 
about 1 month before and after the main spectropolarimet- 
ric runs; a total of 5 and 13 measurements were obtained 
in conjunction with our 2008/2009 and 2007/2008 runs re- 
spectively. Additional pho tometric data fro m the All-Sky 
Automated Survey fASAS. |Poimanskilll997l ) were added to 
improve phase coverage; selecting only grade A exposures, 
and further rejecting statistically deviant points, we are left 
with 12 and 9 supplementary ASAS measurements contem- 
poraneous with our 2008/2009 and 2007/2008 runs. Typical 
rms photometric accuracies are about 50 and 100 mmag for 
CrAO and ASAS data respectively. 



4 VARIATIONS OF PHOTOSPHERIC LINES 
AND ACCRETION PROXIES 

Zeeman signatures with typical peak-to-peak amplitudes of 
about 1% are clearly detected at all times in LSD profiles 
of photospheric lines (see Fig. [T]top right panel); temporal 
variations are also detected but remain moderate, the aver- 
age Stokes V profile over the whole run featuring a roughly 
symmetric shape (with respect to line centre) suggesting the 
presence of a significa nt toroidal field com ponent at the sur- 
face of the star (e.g., iDonati et al.ll2005l ). As a result, the 
line-of-sight projected component of the field averaged over 
the visible stellar hemisphere (called longitudinal field and 
estimated from th e first moment of the Stokes V profile, 
iDonati et allll997l ') is weak, varying from -230 G to -|-70 G 
during the 2008/2009 observing run depending on the epoch 
(see Fig. [2] lower left panel). 

LSD Stokes I profiles of photospheric lines also vary 
with time, both in position and strength (see Fig. [T] top left 
panel). The corresponding RV variations reach a peak-to- 
peak amplitude of about 2 kms~^ (about a mean of about 
17 kms~^), similar to previous results reporting that these 
variations are apparently stable on a long-term basis and 
tentatively attributed to spots at the surface of AA Tau 
l|Bouvier et al.ll2007a ) . By comparing our spectra of A A Tau 
to those of a spectroscopic template of similar spectral type, 
we can in principle retrieve estimates of the veiling, i.e., 
the amount by which photospheric lines are weakened with 
respect to those of the template star. The template star we 
observed (61 Cyg B, of similar temperature but different 
luminosity class) shows LSD profiles that are weaker (by 
10 — 30%) than most of our AA Tau profiles. While it first 
indicates that a better template star is needed to obtain 
absolute veiling estimates, it at least suggests that veiling is 
only moderate for AA Tau at the time of our observations. 



Using our strongest LSD profile of AA Tau (at cycle 53.103) 
as a zero- veiling reference, we find that veiling is most of the 
time smaller than 5% during the 2008/2009 observing run, 
except between rotation cycles 49.7 — 50.2 where it peaks at 
about 20% (see Fig. [2] lower right panel). 

We find that LSD Stokes I and V profiles collected at 
similar rotation phases can exhibit different shapes at differ- 
ent rotation cycles (e.g., at cycles 48.991, 49.961 and 52.989), 
indicating that intrinsic variability is significant. In particu- 
lar, veiling does not repeat well between successive rotation 
cycles (see Fig. [5] lower right panel), suggesting that this 
intrinsic variability likely refiects unsteady accretion at the 
surface of AA Tau. We note that veiling remains smaller 
than 5% throughout our 2009 January observations (ro- 
tation cycle 52.3 — 53.3), suggesting that we are catching 
AA Tau in a state of constantly very low accretion at this 
time of the 2008/2009 run; we suspect that, as a result, this 
is also an epoch at which rotational modulation of spec- 
tral lines is easiest to detect against intrinsic variability. We 
finally note that strong veiling (when present) occurs at ro- 
tation phases 0.8 - 1.1 in the 2008/2009 run, i.e., slightly 
before eclipse maximum (at phase 1.05, see Fig. [2] upper 
panel) and similar to previous reports (e.g., iBouvier et al.l 
l2007bl ). 

Usually considered as the most reliable accretion proxy. 
He I Da emission at 587.562 nm is thought to be produced 
in the postshock region at the footpoints of accretion fun- 
nels, i.e., in chromospheric accretion spots. At the time of 
our observations. He I emission always shows up as a nar- 
row profile (see Fig. [3l top panel) and amounts to an average 
equivalent width of about 15 kms"'^ (i-e., 0.030 nm). It is 
strongly variable with time, with equivalent widths varying 
from 7.5 to 45 kms"^ (i.e., 0.015 to 0.090 nm) throughout 
our 2008/2009 run (see Fig. [2] second panel right column). 
Previous studies (e.g., JBouvicr ct aL 2003, 2007b) report 
very similar levels of He I emission, equal to about 25 km s"'^ 
in average (i.e., 0.05 nm), most of the time lower than 
50 kms~^ (i.e., 0.10 nm) and reaching at most 75 kms"^ 
(i.e., 0.15 nm) during highest-accretion episodes. Times of 
strongest He I emission clearly coincide with epochs of max- 
imum veiling, as pre viously reported in the literature (e.g., 
iBouvier et al.ll2007bl ); in particular, using the published cor- 
relation between vei ling and He l emission (see Fig. 5 of 
IBouvier et alj|2007bf) , we can safely confirm that veiling is 
indeed very small at rotational cycle 53.103 (assumed as 
our zero- veiling reference epoch, see above). Most of the ob- 
served variation of He I emission does not repeat from one 
rotation cycle to the next and thus mostly reflects intrin- 
sic variability rather than rotational modulation. We find 
that He I emission is minimal throughout our 2009 January 
observations, with equivalent widths varying from 7.5 to 
10.2 kms"'^ (i-e., 0.015-0.020 nm); it confirms that rotation 
cycle 52.3 — 53.3 corresponds to a very-low-accretion stage 
of A A Tau with limited (i.e., ~30% peak-to-peak) though 
definite rotational modulation, and with maximum emission 
occurring around phase 1.0. 

Clear Zeeman signatures are detected in conjunc- 
tion with He I Ds emission, corresponding to longitudinal 
fields as strong as —2.5 kC (see Fig. [2] second panel left 
column) and similar t o pre vious published reports (e.g., 
IValenti fc Johns- Krul] \2004 ). We note that longitudinal 
fields are not markedly different in phases of low and high 
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accretion, indirectly confirming that He I emission is mostly 
probing accretion spots at the footpoints of accretion fun- 
nels; as a result, intrinsic variability on this longitudinal 
field estimate is limited while rotational modulation is fairly 
clear. Longitudinal field is strongest around phase 1.0 (i.e., 
in phase with line emission), indicating that this is when the 
accretion spot is best visible to an Earth-based observer. 

Core emission in Ca ll IRT lines are also useful for prob- 
ing magnetospheric accretion. Since Ca ll emission is pre- 
sumably also coming from the non-accreting chromosphere, 
it is less specifically related to accretion spots than He I 
emission and is thus a more ambiguous proxy. However, 
the redder spectral location, higher magnetic sensitivity and 
multiple nature of the corresponding spectral lines more 
than compensate for this drawback; moreover, the shape of 
the corresponding Zeeman signatures is simpler (i.e., fea- 
turing a classical, nearly-antisymmetric pattern, as opposed 
to the strongly non-antisymmetric Zeeman signatures of the 
He I line, e.g., iDonati et al.l |2007| . [20083) ^nd thus easier 
to model. To extract the core Ca ll emission profiles (used 
later in the imaging analysis, see Sec. [5} from our spectra, 
we start by constructing a LSD-like weighted average of the 
3 IRT lines; we then subtract the underlying (much wider) 
photospheric absorption profile, with a Lorentzian fit to the 
far wings over a velocity interval of ±200 kms~^ about the 
emission core. Equivalent widths of the Ca ll emission are in 
average equal to about 15 kms~^ (i-e., 0.043 nm) and vary 
from 12 to 24 kms"^ (i.e., 0.035 - 0.070 nm) during our 
2008/2009 run. As for He I, most of the variation does not 
repeat between successive rotation cycles (see Fig. [T] bottom 
left panel and Fig. [2] third panel right column) and is thus 
intrinsic in nature. Note however that the strongest Ca ll 
emission episode (at cycle 48.510) coincides with no more 
than a small increase in He I emission (see Fig. [2] right col- 
umn) while large He I emission episodes generate only mod- 
erate Ca II emission; this directly relates to the fact that 
Ca II emission is more sensitive to chromospheric activity 
and flares (the likely cause of the emission episode at cycle 
48.510) than He I emission. During rotation cycle 52.3 — 53.3 
(i.e., in 2009 January), Ca ll emission is close to minimum, 
varying by about 20% peak to peak; maximum emission is 
reached at about phase 1.0 (as for He i emission, see Fig. [5] 
middle panels right column), confirming that most of the 
observed rotational modulation at this time of the run is 
due to the chromospheric accretion spot. 

Zeeman signatures in Ca ll IRT lines are clearly de- 
tected at almost all epochs, with average peak-to-peak am- 
plitudes of about 10% (see Fig. [T] bottom right panel). Cor- 
responding longitudinal fields reach up to about —1 kC in 
high accretion states (—0.9 and —1.2 kC at cycles 49.961 
and 50.082 respectively) and more typically range from to 
—700 G otherwise. In 2009 January in particular (rotation 
cycles 52.3 — 53.3), longitudinal fields vary very smoothly, 
and in phase with longitudinal fields derived from the He I 
line; this is further evidence that the observed fluctuation 
mainly reflects rotational modulation. Maximum fleld is 
reached at phase 0.95, i.e., when the accretion spot faces the 
observer; maximum fleld strength is about 3 times weaker 
than that from He I lines, suggesting that the accretion spot 
contributes no more than about one fourth to Ca II emission 
(the remaining part coming from the non-accreting chromo- 
sphere). 



Ha and H/3 both feature emission components of vary- 
ing strengths (see Fig. [3] bottom panel); while Ha emission 
is strong and wide and truncated by a central absorption 
component (giving the profile a mostly double-peak appear- 
ance), H/? emission is much weaker and mostly conflned 
to the blue proflle wing. During the 2008/2009 run, the 
average equivalent width of Ha is 550 kms~^ (1.2 nm), 
with peaks of up to 2050 kms"'^ (4.5 nm, at cycle 48.991). 
For H/3, the average equivalent width is close to (within 
about 20 kms^^ or 0.03 nm), with the weak blue- wing emis- 
sion barely compensating for the central absorption; at cy- 
cle 48.991, H/3 emission peaks at an equivalent width of 
430 kms"^ (0.7 nm). 

In 2009 January, Ha is weakest, with an average equiv- 
alent width of 250 kms~^ (0.55 nm); blue and red emis- 
sion components vary more or less in phase with each other, 
peaking around phase 0.35 and minimum at phase 0.85. At 
the same time, H/3 exhibits strong central absorption and re- 
duced blue- wing emission, with an average equivalent width 
as low as —70 kms^^ (—0.1 nm) and proflle variation con- 
centrating mainly in the red wing (at about -f-lOO kms~^); 
maximum and minimum emission occur at ph ase 0.5 and 
1.0 respectively. Following iBouvier et al.l (l2007br ). we specu- 
late that the absorption episodes in the red wing emission of 
Yip trace accretion funnels crossing the line of sight as they 
corotate around the star. 

Using the average equivalents widths of the He I line, 
the Ca II IRT and Ha during our 2008/2009 run, we derive 
logarithmic line fluxes (with respect to the luminosity of the 
Sun Lq) equal to —5.3, —5.1 and —3.6 respectively. This im- 
plies a logarithmic accretion luminosity (with res pect to Lp) 
of —2 .3 ±0.3 using the empirical correlations of IFang et al.l 
(J2009r) and putting strong weight on the He I data, and thus 
a logarithmic mass-accretion rate at the surface of the star 
of —9.2 ± 0.3 (in M0yr~^). Mass accretion rates can also 
be estimated (though more roughly) through the full width 
of H a at 10% height (e.g., iNatta et all |2004| : ICieza et all 
I2OI0I ): in our case, the full widths we determine range from 
350 kms~^ (in 2009 January) up to 480 kms~^ (at cycle 
48.991), with an average of about 420 kms~^, yielding log- 
arithmic mass-accretion rates of —9.5, —8.8 and —8.2 (with 
typical errors of ±0.6) respectively supporting our main es- 
timates. We find that the logarithmic surface mass-accretion 
rate at magnetic maximum (i.e., around phase 0.0) can vary 
by about an order of magnitude (from —9.4 to —8.5) from 
one cycle to the next, while that at magnetic minimum (i.e., 
around phase 0.5) is more stable on the long term (at about 
—9.6). During the low-accretion phase of 2009 January, the 
logarithmic mass-accretion rate at the surface of AA Tau 
is found to be —9.5 in average and to vary by only ~ 50% 
peak-to-peak (from —9.6 to —9.4). 

For our 2007/2008 run, we obtain similar results (not 
shown here) but with larger error bars and worse phase sam- 
pling (7 phases only, the first 6 roughly covering one cycle, 
see Table [2]). Zeeman signatures in Stokes V LSD profiles 
of photospheric lines are again detected at all phases, trac- 
ing longitudinal fields of ±90 to —220 G. Emission in Ca 11 
and He I lines vary smoothly with time, reaching maximum 
strengths at phase 1.3; longitudinal fields from the Ca 11 
lines also show smooth changes from to —730 G, reach- 
ing maximum field at phase 1.2. This suggests in particu- 
lar that the accretion spot is located at phase 1.2-1.3, i.e.. 
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Figure 1. Temporal variations of the Stokes / (left) and Stokes V (right) LSD profiles of the photospheric lines (top) and of the Ca II 
emission (averaged over the 3 IRT lines, bottom) of AA Tau in 2008 December (left column of each panel) and 2009 January (right 
column). The Stokes / profiles of Ca II emission is shown before subtracting the underlying (much wider) photospheric absorption, hence 
the reduced flux in the far wings (with respect to the unit continuum). To emphasize variability, the average profile over the run is shown 
in red. Rotation cycles (as listed in Table 1) and 3a" error bars (for Stokes V data only) are also included next to each profile. 
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Figure 2. Temporal variations of the integrated brightness (top panels), He I Ds (second panels), Ca II IRT (third panels) and photo- 
spheric lines (bottom panels) of AA Tau in 2008 December and 2009 January. Longitudinal field variations are shown on the left side 
(3 lower panels) while equivalent width and veiling variations are shown on the right side. Data collected within cycles 48.3 — 49.3, 
49.3 — 50.3 and 52.3 — 53.3 are respectively plotted as red circles, green squares and blue triangles. Fits with sine/cosine waves (plus 
first overtones for the 2 lower left panels) are included (and shown as dashed lines) for data collected within rotation cycle 52.3 — 53.3 
(to outline the amount of variability caused by rotational modulation) and for all photometric data; ±1 cr error bars on data points are 
also shown whenever larger than symbols. 



again slightly before eclipse maximum (i.e., brightness min- 
imum, occurring at phase 1.35). The variability observed in 
all spectral proxies is apparently compatible with rotational 
modulation; the moderate phase coverage and poor phase 
redundancy achieved in this initial run does however not al- 



low us to quantify reliably the relative strength of intrinsic 
variability. Emission in He I, Ca ll IRT and Ha lines fea- 
ture equivalent widths of 14 kms~^ (0.027 nm), 14 kms~^ 
(0.040 nm) and 700 kms"'^ (1.5 nm), corresponding to an 
average logarithmic mass-accretion rate of —9.3 ± 0.3 (in 
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Figure 3. Variation of the He I D3 (Stokes /: top left, Stokes V: top right), Ho (bottom left) and H/3 (bottom right) lines of AA Tau 
in 2008 December (left column of each panel) and 2009 January (right column). Line fluxes in the far wings of all Stokes / profiles are 
close to the (unit) continuum level. To emphasize variability, the average profile over the run is shown in rod. Rotation cycles (as listed 
in Table 1) and 3cr error bars (for Stokes V data only) are also mentioned next to each profile. 
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Mq yr"^) similar to that during the 2008/2009 run. The log- 
arithmic surface mass-accretion rates at magnetic minimum 
and maximum are respectively equal to —9.6 and —9.0, i.e., 
showing a larger peak-to-peak variation (but a similar level 
at magnetic minimum) than in 2009 January. 



5 MODELLING THE SURFACE OF AA TAU 

The visual inspection of all spectral proxies and their tempo- 
ral variations (presented in SecQ straightforwardly reveals 
a number of interesting results, in particular concerning the 
chromospheric spot at the footpoint of the accretion funnels 
linking AA Tau with its accretion disc. Assuming that the 
field within this accretion spot is mainly radial and gener- 
ates the periodic longitudinal field fluctuation seen in the 
He I line, we can derive from the observed rotational mod- 
ulation that the accretion spot is located at a latitude of 
about 70° (to produce a null field extremum at phase 0.5 
and given that i ~ 70° , see Sec. [2| and thus that the field 
strength within the spot is of order of 3 kG (generating lon- 
gitudinal field maxima of about 2 kG when the field reaches 
its minimum inclination of 50° to the line of sight). Eclipse 
maxima are obviously linked with (though slightly lagging) 
longitudinal field maxima and peak absorption in the red 
wing of Balmcr lines, confirming that eclipses are likely due 
to a disc warp at the base of the accretion veil, as initially 
pr oposed bylBouvier et al.l (|2003l ) and further documented 
bv lBouvier et al.l l|2007bl ). 

Going further, and in particular quantifying how mag- 
netic fields are distributed at the surface of AA Tau and 
investigating how fields connect the protostar to its accre- 
tion disc, requires a more detailed modelling of the observed 
profiles and their rotational modulation. This is what we 
propose below where we introduce a more mature and elab- 
orate model that builds up on the prelimi nary version intro- 
duced in previous papers (jDonati et al.l 12007 . 2008al ) . We 
then apply the model to our spectra of AA Tau collected 
in 2009 January (those from 2008 December containing too 
much intrinsic variability, see Sec. [3| and in 2007/2008. 



5.1 Mapping magnetic field, photospiieric 
brigiitness and ciironiospheric accretion 

Our model is designed to recover simultaneously maps of 
the surface magnetic field, of the photospheric brightness 
and of the accretion-powered excess emission from sets of 
Stokes 7 and V LSD and Ca ll profiles of classical T Tauri 
stars such as those shown in Fig. [l] We do not attempt at 
matching He I profiles in this paper, as the complex shape 
of their Zeeman signatures (see Sec. |3| would require some 
modelling of the velocity flows within the postshock region of 
accretion funnels, beyond the scope of our study. Similarly, 
a quantitative modelling of Balmer line profiles, requiring a 
full 3D description of the magnetosphere (including density 
and velocity fields), is postponed for future papers. 

In this model, the magnetic field is described through 
its poloidal and toroidal co mponents expressed as spherical- 
harmonics (SH) expansions l|Donati et al.ll2006l ). The spatial 
distribution of photospheric brightness (with respect to the 
quiet photosphere) and that of Ca ii excess emission (with 
respect to the quiet chromosphere) are modelled as series 



of independent pixels (typically a few thousand) on a grid 
covering the visible surface of the star; brightness spots are 
assumed to be darker than the photosphere while accretion 
spots are supposed to be brighter than the chromosphere. 
The main difference with respect to the preliminary version 
of the code is that the photospheric and chromospheric im- 
ages are no longer assumed to be homothetic but are fitted 
to the data independently from one another; moreover, the 
quantity that we now recover for the photospheric image is 
brightness (rather than spottedness) . 

Following the principles of maximum entropy, the code 
automatically retrieves the simplest magnetic topology, 
brightness image and accretion map compatible with the 
series of rotationally modulated Stokes / and V LSD pro- 
files. The reconstruction process is iterative and proceeds 
by comparing at each step the synthetic Stokes I and V 
profiles corresponding to the current images with those of 
the observed data set. To compute the synthetic profiles, 
we add up the elementary spectral contributions from each 
image pixel over the visible hemisphere, taking into account 
all relevant local parameters of the corresponding grid cell 
(e.g., brightness, chromospheric excess emission, magnetic 
field strength and orientation, radial velocity, limb angle, 
projected area). Since the problem is partly ill-posed, we 
stabilise the inversion process by using an entropy criterion 
(applied to the SH coefficients and to the brightness/excess 
emission image pixels) aimed at selecting the field topolo- 
gies and images with minimum information among all those 
compatible with the data. The algorithm used to solve the 
maximum entropy problem has be en adapted from that pre- 
sented bv lSkilling fc BrvanI ([1984') which iteratively adjusts 
the image by a multidirection search in the image space; in 
particular, this algorithm was found to be much more effi- 
cient than standar d conjugate gradient techniques (see, e.g., 
iBrown et al.lll99ll . for more details). The relative weights 
attributed to the various SH modes can be imposed, e.g., 
for purposedly producing antisymmetric or symmetric field 
topologies with respect to t he centre of the star (by f avour- 
ing odd or even SH modes. 'Donat i et alllJOOTl , 1200831 ). 

The local synthetic photospheric Stokes I and V line 
profiles emerging from a given grid cell (noted Jp and Vp) 
are modelled using the following equations: 



^P = 



6 (V'J'p,. 



+ (1 - V')/p,q) 



(2) 



where /p,ni and Jp^q are the local Stokes / photospheric pro- 
files corresponding to the magnetic and non-magnetic areas, 
14>,m the local Stokes V photospheric profile corresponding 
to the magnetic areas, fe (0 < 6 < 1) the local brightness 
relative to the quiet photosphere and V' (0 < V' ^ 1) the 
relative proportion of magnetic areas within the grid cell 
(called filling factor) . For simplicity, we further assume that 
/p,m and 7p,q differ by no more than magnetic effects, i.e., 
that 7p,q equals /p,m taken at B = 0. 

Similarly, we describe the Stokes / and V Ca ll emission 
profiles (noted /o and Vc) with equations; 

h = /k + e/a = V'(-^k,m + e/a,m) + (1 - V')(/k,q + e/a,q) 
Ve = li -H eVa = l/'(T4,m + eVa,m) (3) 

where /k and Vk on the one hand, and Ja and 14 on the 
other hand, are the respective contributions to / and V pro- 
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files from the quiet chromosphere and the accretion regions, 
and e (e > 0) the local excess emission from accretion re- 
gions or equivalently the fraction of the grid cell occupied by 
accretion regions; the contributions from magnetic regions 
to /k and /a (respectively 14 and 14) are denoted Jk,m and 
Ia.,m (respectively 14, m and Vk.m) while the contributions 
from non-magnetic regions are denoted /k,q and 7a, q. For 
simplicity, we further assume that 7a and 7k (and similarly, 
7a,m and 7k,m, 7a,q aud 7k,q, 14 and 14, 14,m and 14,m) are 
homothetic, with a scaling factor denoted e (e > 1); we also 
suppose that 7k, m and 7k, q (and thus 7a,m and 7a, q) differ 
by no more than magnetic effects (just as 7p,ni and 7p,q). 
The filling factor tp is finally assumed to be the same for 
the photospheric, chromospheric and accretion profiles, and 
for all grid cells. Local profiles from all grid cells are then 
integrated taking into account the cell visibility, projected 
area, limb-angle and radial velocity. 

To describe the local profiles (i.e., 7p,m, 7p,q, 14,m, 
7k,m, 7k, q, 14, m), we use Unno-Rachkovsky's equations 
known to provide a good description of Stokes 7 and 
V profiles (including magneto-optical effects) in the pres- 
ence of both weak and strong m agnetic fields (e.g., 
iLandi degl'Innocenti fc Landolfill2004l . Sec. 9.8) despite their 
being based on the assumption of a simple Milne-Eddington 
atmospheric model. For photospheric LSD profiles (i.e., 7p,ni, 



14. 



we set the central wavelength, Doppler width 



and Lande factor of our equivalent line to 640 nm, 2 kms^^ 
and 1.2 respectively and adjust the average line equivalent 
width to the observed value (with veiling removed from all 
profiles); for Can chromospheric profiles (i.e., 7k, m, 7k. q, 
14, m), the central wavelength, Doppler width and Lande 
factor are set to 850 nm, 6.7 kms~^ and 1 respectively, the 
equivalent width being set to slightly below the minimum 
Ca II emission in our observations. 

By comparing synthetic and observed profiles, we derive 
the spatial distributions of the local photospheric brightness 
b, of the Ca ll excess emission e due to accretion spots, and 
of the magnetic vector B over the stellar surface. As a side 
product, we also obtain estimates of several other model 
parameters (by selecting the image with minimum informa- 
tion content at a given fit accuracy), and in particular of 
the line rotational broadening vsini, of the RV Wrad, of the 
filling factor ^p and of the accretion profile scaling factor e. 
Simulations demonstrate that the code can reliably recon- 
struct magnetic topologies, brightness images and accretion 
maps; in particular, it is found to be very efficient at recov- 
ering the poloidal and toroidal field components, a very use- 
ful diagnostic when studying large-scale magnetic topologies 
produced by dynamo processes at the surface of cool stars. 
Surface differential rotation patterns shearing the bright- 
ness and magnetic images can also be reliably retrieved, as 
demo nstrated in its most recent application (JDonati et al.l 
|2010D . 

5.2 Application to A A Tau 

The imaging model described above assumes that the ob- 
served profile variations are mainly due to rotational mod- 
ulation, and possibly to surface differential rotation as well 
when the star is observed for at least several rotation cy- 
cles; all other sources of profile variability (and in par- 
ticular intrinsic variability like flaring or short-term high- 



accretion episodes) cannot be properly reproduced and thus 
contribute as noise into the modelling process, degrading 
the imaging performance and potentially even drowning all 
relevant information. 

Filtering out significant intrinsic variability from the 
observed profiles of AA Tau is thus worthwhile to optimise 
the behaviour and convergence of the imaging code, lead- 
ing us to discard our 2008 December data and to concen- 
trate only on those collected in 2009 January and during 
our 2007/2008 run. For the latter 2 sets, we further need 
to suppress veiling, e.g., by scaling all LSD Stokes 7 and V 
photospheric profiles to ensure that unpolarized lines have 
the same equivalent widths; we also need to retain rotational 
modulation only in the variations of Ca ll IRT profiles, e.g., 
by fitting them with a sine-|-cosine wave and scaling the 
corresponding Stokes 7 and V profiles to ensure that unpo- 
larized lines match the fitted equivalent widths (see Fig. |4]). 
While obviously no more than approximate, this procedure 
has the advantage of being very straightforward yet reason- 
ably efficient , and proved succes sful in the case of the cTTS 
V2247 Oph (JDonati et al.ll2010l '). For both runs, rotational 
modulation in Ca ll emission is found to be of order 30% 
peak-to-peak. 

Given the moderate v sin i of AA Tau, we limit the SH 
expansions describing the field at ^ = 5. The reconstructed 
magnetic, brightness and accretion map^J of A A Tau are 
shown in Fig. [5] for both epochs, with corresponding fits 
to the data shown in Fig. [6] The overall fits are good, re- 
producing the data down to the noise level starting from 
initial reduced chi-squares of 18 and 5.5 for the 2009 Jan- 
uary and 2007/2008 data sets respectively. Observations at 
both epochs (covering roughly a single rotation cycle) do 
not allow an estimate of surface differential rotation. Op- 
timal fits are obtained for vsini = 11.5 ± 0.5 kms~^ and 
Wrad = 17.2 ± 0.1 kms~^, in good agr e ement with previ- 
ous estimates (e.g., iBouvier et aLll2003l . l2007bh : Ca ll IRT 
lines are found to be slightly redshifted (by about 1 kms~^) 
with respect to photosphe ric lines, as usual for cTTSs (e.g., 
iDonati et alll2007l . [JOOSah . We also find that filling factors 
of about i/) ~ 0.3 ±0.1 provide the best fit to the far wings 
of the observed Stokes V profiles (see Fig. |6]). 

The emission profile scaling factor e, describing the 
emission enhancement of accretion regions over the quiet 
chromosphere, is set to e = 10. This value is however some- 
what arbitrary; fitting the observed strengths of emission 
profiles fixes e'^e rather than e alone, with e describing 
how the excess emission from accretion regions varies over 
the surface of the star (see Eq[3| and ^ e thus estimating the 
fractional area of accretion footprints. We nevertheless think 
that our choice of e is realistic (and accurate within a typical 
factor of a few) given that the fractional areas of accretion 
footprints it yields (a few %, see belo w) are grossly compat- 
ible with published estimates (e.g., IValenti fc Johns-Krulll 
|2004| ). 

The reconstructed large-scale magnetic topologies are 



^ Reconstructed maps are only shown down to latitudes of —30°; 
being both strongly limb-darkened and visible for only a short 
amount of time, features at lower latitudes contribute little to 
the observed profiles and are therefore basically out of reach of 
imaging techniques, especially when phase coverage is moderate. 
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Figure 4. Measured (red open circles) and fitted (blue pluses) equivalent widths of the Ca II IRT LSD profiles throughout our 2009 
January (left) and 2007/2008 (right) observing runs. The model wave (green line) providing the best (sine+cosine) fit to the data (with a 
period of 8.22 d) presumably traces rotational modulation, while the deviation from the fit illustrates the strength of intrinsic variability. 



similar at both epochs and essentially consist in a ~ 2 kG 
dipole inclined at about 20° to the rotation axis, concen- 
trating about 90 — 95% of the poloidal field energy. This is 
about 30% weaker than predicted from the He i alone (indi- 
cating a polar field of about 3 kG, see above), likely refiect- 
ing the uncertainty in the modelling of the Ca ll profiles. We 
thus conservatively conclude that the dipole component of 
A A Tau is in the range 2 — 3 kG. These estimates are grossly 
compatible with previou s results from dipolar fits to He I 
spectropolarimetric data (|Valenti fc Johns-Krullll2004l ). The 
quadrupole/octupole components are significantly weaker, 
by typically a factor of 5 to 10 depending on the relative 
weights given to the different SH modes (see Sec.[5ll. We also 
find that the field includes a significant toroidal component 
in the shape of an equatorial ring of negative (i.e., clock- 
wise) azimuthal field (0.5 — 1 kG), totalling about 15 — 20% 
of the magnetic energy at the surface of the star. As al- 
ready mentioned in Sec. |31 this toroidal component can be 
straightforwardly traced back to the average shape of LSD 
Stokes V profiles of photospheric lines (mostly symmetric 
about the line centre, see Fig.[T| and can thus be considered 
as reliable. 

The main difference between the two reconstructed 
magnetic images is an apparent shift of about 0.25 cycle be- 
tween the phases towards which the dipoles are tilted (phase 
in 2009 January and phase 0.25 during the 2007/2008 run). 
Since eclipse phases (respectively equal to 0.05 and 0.35, see 
Sec. 131) are also shifted by about the same amount, we can 
conclude that this shift essentially refiects the uncertainty 
on the rotation period; a slightly smaller period (of ~ 8.18 d 
instead of 8.22 d, s till compatible with the quoted uncer- 
tainty of 0.03 d of iBouvier et"aLl l2007bh would bring the 
two images (and photometric data sets) in phase with one 
another. Given the moderate phase coverage and the level 
of intrinsic variability observed on AA Tau (even in Stokes 
V profiles of photospheric lines, see Fig. (21), it is difficult 
to ascertain whether the remaining differences in the recon- 
structed large-scale fields at both epochs are real. 

The surface brightness maps we reconstruct are also 
grossly similar once accounting for the phase shift between 
both epochs. They mostly feature a large spotted area cen- 
tred at latitude ~ 50° , covering about 10% of the full surface 



and spreading over about half the full phase range. This dark 
area accounts in particular for the observed Stokes I profile 
asymmetries (see Fig. [§] left panels) and for the correspond- 
ing RV variations. The location of this spot agrees well with 
the epochs of eclipse maxima, as expected from RV max- 
ima and minima occurring respectively ~0.25 rotation cycle 
before and after eclipse maxima (see Fig. (Tjl . 

Given the depth of eclipses (ranging from 0.5 to 
1.5 mag), it is fairly clear that this dark region can- 
not fully account by itself for the observed photometric 
variability; color and photopolarimetric variations further 
demonstrate that circumstellar extinct ion is the main cause 
of the observed b rightness changes (jBouvier et al.l l2003l : 
iMenard et al.ll2003l ). Similarly, one can wonder whether the 
reconstructed dark region is not spurious, i.e., whether the 
observed line profile variations are not the direct conse- 
quence of eclipses. This is unlikely for at least two reasons. 
Firstly, the varying asymmetries in the observed Stokes / 
photospheric profiles and the corresponding RV fiuctuations 
are known to occur at all time s, even when eclipses are miss- 
ing (e.g., around JD 2.451.522. iBouvier et al.ll2003l ): our own 
data also demonstrate that RV variations are grossly stable 
and weakly affected by changes in the accretion rate and veil- 
ing strength (see Fig. [7] and bottom right panel of Fig. [5|. 
Secondly, similar (though smaller) dark regions, more or 
less spatially coincident with m agnetic poles, are ob served 
on similar cTTSs (e.g., BP Tau. iDonati et al.ll2008"al ). even 
when not eclipsed by a warped inner disc (as AA Tau). We 
thus conclude that the cool region reconstructed at the sur- 
face of AA Tau is most likely real, at least its darkest parts 
found to grossly overlap with magnetic poles; it is however 
possible that some of it (e.g., the trailing extension towards 
phase 0.45 in the 2007/2008 image, see lower panels of Fig.[5| 
is spurious and caused by minor line profile changes directly 
caused by eclipse episodes. 

The accretion maps we derive at both epochs include 
one high-latitude region covering about 2% of the total stel- 
lar surface (when assuming e = 10, see above) and shifted 
again by about 0.25 cycle between both runs. This accre- 
tion region is crescent-shaped, i.e., elongated along paral- 
lels and spanning a phase range of about 0.2. The brightest 
spot within this accretion region is leading, crossing the line- 
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Figure 5. Maps of the radial, azimuthal and meridional components of the magnetic field B (first and third rows, left to right panels 
respectively), photospheric brightness b and excess Ca II IRT emission e (second and fourth rows, first and second panels respectively) 
at the surface of AA Tau, during our 2009 January (two upper rows) and 2007/2008 (two lower rows) runs. Magnetic fluxes are labelled 
in G. In all panels, the star is shown in flattened polar projection down to latitudes of —30°, with the equator depicted as a bold circle 
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Figure 6. Maximum-entropy fit (tliin red line) to the observed (thick black line) Stokes / and Stokes V LSD photospheric profiles (first 
two panels) and Ca II IRT profiles (last two panels) of AA Tau, for our 2009 January (upper row) and 2007/2008 (lower row) data 
sets. The light-blue curve in the third panels shows the (constant) contribution of the quiet chromosphere to the Stokes / Ca II profiles. 
Rotation cycles and 3(T error bars (for Stokes V profiles) are also shown next to each profile. 
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Figure 7. RV variations of phiotospheric lines of AA Tau in 2008 
December and 2009 January. Data collected within cycles 48.3 — 
49.3, 49.3 — 50.3 and 52.3 — 53.3 are respectively plotted as red 
circles, green squares and blue triangles, with a sine/cosine fit 
to the blue triangles shown (dashed line) to outline rotational 
modulation. 



of-sight slightly before epochs of eclipse maxima (by about 
0.1 cycle); it is also roughly coincident with the visible mag- 
netic pole. These crescent-shaped accretion spots are very 
similar to those predicted in numerical simulation of mag- 
netos pheric accretion to sligh tly inclined magnetic dipoles 



(e.g.. iRomanova et al]l2004bl) . 



6 SUMMARY AND DISCUSSION 

This paper presents the first results of the MaPP project 
aimed at studying the impact of magnetic fields on star and 
planet formation. In particular, we report the detection of 
Zeeman signatures from the prototypical cTTS AA Tau at 
two different epochs (2008 December / 2009 January and 
2007 December / 2008 January) and using the ESPaDOnS 
and NARVAL spectropolarimeters on CFHT and TBL re- 
spectively. From phase-resolved data sets, we successfully 
derive maps of the magnetic field, of the surface bright- 
ness and of the accretion-powered emission on AA Tau. We 
briefly summarize below the main results and discuss their 
implications for our understanding of magnetized stellar for- 
mation. 

The large-scale magnetic field we reconstruct consists 
mostly of a 2 — 3 kG dipole slightly tilted (by about 20°) 
to the rotation axis, with higher SH modes of the poloidal 
component being at least 5x smaller. This is different from 
BP Tau, the only other cTTS of similar mass magnetically 
imaged to date, where the dipole and octupole poloidal com - 
ponents had comparable strengths (jDonati et al.l l2008al ). 
This is even more radically different from those of both 
lower- and higher-mass cTTSs (namely V2247 Oph and 
MT Ori) for which the la rge-scale poloidal f ield is much 
weaker and more complex IjDonati et al.ll20ld . Skelly et al, 
2010, in preparation); similar conclusions we re reached for 
the h igh-mass cTTSs CV Cha and CR Cha (jHussain et al.l 
I2OO9I ). This makes AA Tau grossly similar to mid-M dwarfs 
(whose large-scale fields ar e also found to be strong, mostly 
poloidal and axisvmmetric. lMorin et al.ll2008l 'l. while higher- 
and lower-mass cTTSs are closer to early- and late-M dwarfs 



respectively (featuring weak and complex large scal e fields 
in average, iDonati et all l2008d : iMorin et aiTl2010l '). This 
strengthens the idea that magnetic fields of cTTSs are pro- 
duced through dynamo pro cesses like those of M dwarfs 
(JDonati fc Landstreetll2009l ). The origin of the differences 
between the large-scale topologies of AA Tau and BP Tau 
(showing very different dipole to octupole intensity ratios, 
> 5 and ~1 respectively) are unclear, both stars featuring 
very similar masses, ages and rotation rates; we suggest that 
they may reflect different magnetic states in a long-term 
magnetic cycle. 

AA Tau also includes a small but significant toroidal 
component totalling about 15 — 20% of the reconstructed 
magnetic energy, i.e., about twice as much as that of BP Tau. 
This is however different from mid-M dwarfs for which the 
toroidal compon ent rarely exceeds a few % of the overall 
magnetic energy (jMorin et al.ll2008t ). We speculate that this 
toroidal component relates to the accretion flow, although 
the link between both is not clear yet; one option is that 
accretion is slowing down the surface of the cTTS through 
star/disc magnetic coupling (see below), producing a sub- 
surface shear, an additional interface dynamo (on top of the 
underlying dominant distributed dynamo operating within 
the bulk of the star) and the associated toroidal component. 
Regular magnetic monitoring of cTTSs is required to inves- 
tigate this idea further. 

Using the He I D3, Ca 11 IRT and Ha emission fluxes 
(with stronger weight on the flrs t and more reliabl e proxy) 
and the empirical correlation of iFang et al.l (|2009l ) . we de- 
rive that the average logarithmic accretion rate at the sur- 
face of AA Tau (in M© yr"'^) is —9.2 ± 0.3; at magnetic 
maximum (i.e., phase 0.0), the logarithmic surface mass- 
accretion rate is found to vary by about an order of magni- 
tude (from —9.4 up to —8.5) between rotation cycles, while 
it is roughly constant (at about —9.6) at magnetic minimum 
(i.e., phase 0.5). Given the estimated strength of the dipole 
held of AA Tau (2 - 3 kG at the pole, i.e. 1.0 - 1.5 kG at 
the equator, see above) and assuming an average logarith- 
mic mass-accretion rate of —9.2, we obtain that the radius 
'"mag at which the inner disc of AA Tau is truncated should 
be equal to rmag — 15 Ri, or equiv alently that r^^^jr^p -c ~ 2 
(using the theoretical estimates of lBessolaz et al.ll2008l . with 
B^ ~ 1.2 kG). 

This result is in conflict with the conditions under 
which accretion can proceed, requiring rmag < fcor to en- 
sure that the effective gravity (i.e., the sum of gravitational 
and centrifugal forces) along accretion funnels is pointing 



towards the star; r-n 



> Tc. 



implies indeed that effec- 



tive gravity points outwards in the inner disc regions, and 
thus that the magnetic star is in the propeller regime (e.g., 
iRomanova et al.l l2004al ) and cannot accrete material from 
the disc, in contradiction with observations. In practice, the 
situation is likely more complex than this simplistic picture; 
numerical simulations of star/disc magnetic coupling in the 
propeller regime (and with rmag ^^ rcor) show for instance 
that a small fraction of the disc material can succeed at 
finding its way into the closed magnetosphere and onto the 
stella r surface even though most of it is expelled o utwards 
(e.g.. IRomanova et al.ll2005l : IUstvugova et ahlbOOd ). In this 
context, the mass-accretion rate at the surface of the star 
(i.e., the one that we estimate from emission proxies) is only 
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a lower limit to the mass-accretion rate in the inner disc re- 
gions. 

The logarithmic accretion rate in the disc needed to 
enforce rmag ^ rcor is ranging from —8.3 to —8.0 (depend- 
ing on the assumed value of B*), i.e., 1.6 — 3.2 times larger 
than the maximum accretion rate observed at the surface 
of the star and 8 — 16 times larger than the average es- 
timated surface accretion rate. Interestingly enough, loga- 
rithmic accretion rates at the surface of the star derived 
from the strongest He I emission levels typically reported 
in t he literature (equal to abou t 75 kms~^ or 0.15 nm, 
e.g., iBouvier et al.l |2003| . l2007bl ) and tracing the highest- 
accretion phases of AA Tau, are equal to —8.2, i.e., compat- 
ible with the mass-accretion rate in the disc that we inde- 
pendently derive from matching rmag and rcor. It confirms 
at least that the upper limit in the surface accretion rate 
that our study predicts is grossly compatible with published 
observations; it also suggests that A A Tau is in a state where 
most of the material accreting through the disc is expelled 
outwards (the propeller regime) and only a small amount (up 
to 40% in 2008 December, ~ 10% in average, ~ 5% in 2009 
January) accreted towards the star. Our observations are in 
this respect qualitatively similar to results of numerical sim- 
ulations of cTTSs in propell er regime (e.g., lRomanova et al.l 
I2OO5I : lUstvugova erahl boogl . 

The strong variability that mass accretion at the sur- 
face of AA Tau is subject to is another argument favouring 
our interpretation. It suggests in particular that the accre- 
tion variability observed for AA Tau mostly relates to the 
variable efficiency at which the disc material succeeds at en- 
tering the closed magnetosphere in propeller regime, rather 
than to an intrinsic variability of the accretion rate within 
the inner disc. This is also in qualitative agreement with 
results of numerical modelling of the propeller regime (e.g., 
[Romanova ct al. 2005) where the accretion fiow from the in- 
ner disc to the surface of the star never reaches steady state 
and remains intrinsically variable. 

The systematic lag of photometric eclipses with respect 
to magnetic/accretion maxima (see Sec. |3| is another inde- 
pendent element favouring our schematic model. Previous 
observations of AA Tau with the optical and UV monitors 
of XMM-Newton show a behaviour similar to that reported 
here, with UV maximum (presumably coinciding with mag- 
netic and a ccretion max ima) occurring well before the opti- 
cal eclipse (JGrosso et al. 2007 ) , suggesting that this is a reg- 
ular phenomenon. We propose that this time delay between 
magnetic poles and eclipse times indicates that field lines 
connecting AA Tau to its inner accretion disc are twisted as 
a result of their different rotation rates, with leading funnel 
footpoints at the surface of AA Tau and a trailing accretion 
warp (producing the lagging variable eclipse) at rmag ^ rcor. 

Our study brings fresh evidence that, thanks to its 
strong dipole field component, AA Tau is still mainly in 
propeller regime, with most of the material in the inner re- 
gions of the accretion disc being expelled outwards and only 
a small fraction accreted towards the star. It suggests in 
particular that spinning down of cTTSs through star/disc 
magnetic coupling can potentially still be efficient at ages of 
about 1 — 2 Myr. AA Tau thus appears as an optimal labora- 
tory for studying in more details the spinning down of cTTSs 
in propeller regime and for testing predictions of numerical 
simulations on this issue. With similarly detailed analyses 



on a dozen of cTTSs, MaPP should soon bring important 
new material for our understanding of magnetospheric pro- 
cesses and their impact on the angular momentum evolution 
of forming Suns. 
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